Low-energy Effective Theory for Spin Dynamics of Fluctuating Stripes 
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We derive an effective Hamiitonian for spin dynamics of fluctuating smectic stripes from the t-J 
model in the weak coupling limit t > J. Besides the modulation of spin magnitude, the high energy 
hopping term would induce a low-energy anti-ferromagnetic interaction between two neighboring 
"blocks of spins". Based on the effective Hamiitonian, we applied the linear spin-wave theory and 
found that the spin-wave velocity is almost isotropic for La2- x Sr x Cu04 unless the structural effect 
is considered. The intensity of the second harmonic mode is found to be about 10% to that of the 
fundamental mode. 
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There's still a lot of interest on stripe physics 0, 0] • 
The stripe modulation was found in cuprate supercon- 
ductor (3 of the sample La2- x Sr x Cu04 (LSCO) in the 
low-temperature orthorhombic (LTO) structure. It is 
also consistent with the incommensurate magnetic peaks 
observed in the inelastic neutron scattering experiments 
0, in which the peak shift from (ir, tt) towards (tt, 0) with 
the derivation very near to 2nd at doping concentration 
S. By partial substitution of Nd for La, the LTO lat- 
tice structure is distorted to the low-temperature tetrag- 
onal (LTT) structure Q, in which the horizontal stripe is 
enhanced and the fluctuating (dynamic) stripe becomes 
more ordered (static). 

Theoretically, by the mean-field analysis of the single- 
band Hubbard model 0, there is a possibility that the 
stripe phase is formed. On the other hand.by employ- 
ing the Schwinger-boson mean-field theory pjto the t-J 
model, it is found that the spiral spin state y| can also 
give the deviation from (tt, tt) upon doping. However, it 
was studied by many approaches that the uniform phase 
is unstable towards phase separation in the range of inter- 
est ratio t/J 9J. One of the consequence of the phase sep- 
aration is to form stripes, in order to be consistent with 
the neutron scattering experiments. Because of the stripe 
fluctuation around the hole domain, the spins across the 
hole domain should be anti-parallel. By considering just 
a small transverse fluctuation of the stripe, it was found 
that two neighboring spins across the hole domain feels 
an anti-ferromagnetic interaction of coupling J 

Although there is still controversy about the stability 
of stripes [llj . in this paper, the fact that the fluctu- 
ating smectic stripe phase is stable is our assumption. 
By including the Coulomb repulsion (which is neglected 
in the t-J model), the holes can only phase separate at 
microscopic scale instead of full phase separation. In or- 
der to balance the hole kinetic energy and the Coulomb 
repulsion, stripes is the simplest solution among those 
inhomogeneous states at microscopic scale. 

In the following, we are going to derive the low-energy 
spin dynamics from the 2D t- J model of the site-centered 



stripe along y-direction. The period of hole domain is 2R 
where R 3> 1 in our model. 

The hopping term H t in the t-J model is written as 



H f = Ht + H\ 



Hi 



(1) 



where Hf- is the transverse fluctuation of vertical stripes, 
and h} is the ID kinetic motion along the vertical 
stripes. Now 
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+ 4+Rx, a ,a C i+(R-iy£, a ,a) +h.C. (2) 



where a label the smectic stripes and i x is summed over 
the super-cell of period 2R. And also 



(3) 



In the weak coupling limit, t J, the fluctuating 
stripe induces a charge-density wave along a;-direction. 
There's two almost degenerate ground states for Hj*~' a , 
|V>s) and |i/'t), which are expressed as 
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V2 
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where 



Hi) = ao|--4T°4T->+ai|--4T4°T-> 

+a-i\-- | o HI -) + ... 

|V> 2 ) = «o|--T4°T4-)+ai|--T4T°4-) 
+o_i| ••T°4T4 ••> + •■ ■ 
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The basis of \tfji) and \tp2) are different by the up and 
down spins. The corresponding coefficients are equal. 
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Two subspaces {|V'i)} and {^2)} are orthogonal to each 
other. By Eq.Q, which is an orthogonal transformation 
from \tpi) and {ifa) to \ip$) and IV't), two subspaces Hs 
and Ht spanned by {IV's)} and {IV't}} respectively are 
also orthogonal to each other. Eq.© can be written in 
the following matrix form under the above basis, i.e., 



Ht a = (-t) 



/0110 
10 1 
10 
10 

V: i i 



)(-*) 



/0110 

10 1 
10 
10 
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where the two identical (2R + l)-dimensional matri- 
ces correspond to two orthogonal subspaces spanned by 
{IV's)} and {IV't)}- Simultaneously, the Heisenberg term 
in the t-J model is written in the matrix form 



Hi 



<-¥) 
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10 
10 
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)(t) 



/o 
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(1) 

s 

(1) 



(^s|tfj'>s) 



0(J 2 /t) 

-o(J 2 /t) 



(10) 

(11) 



and no correction to the wavefunction since two sub- 
spaces are orthogonal. Solving the eigenproblem for 



P (x) = |(^ s )| 2 = |Wt)| 2 

= tt2a( 1+C0S W) 



(12) 



p{r) = 5 + pi cos(2(f • r) + Pi cos(4<f • f) + 



(13) 



4rt 



are split by an energy difference E T 
the energy difference between two "block of spins" is 



J' = 2RS(E { T 1) - E^>) 



2J5R 2 
R+l 



(14) 



Substitution of q = 2ir6 gives J' = J/(2(l + 46)). No- 
tice also that the energy difference per site between anti- 
phase and in-phase becomes J'/(2R) = JS/(1 + AS), 
which is almost linear in 5 for S -C 1. For S = 1/8, energy 
difference is J/12 ~ 130K (take J = 135meV), which is 
consistent with the observation of the incommensurate 
peak up to around 100-fT p|. 

There is then an effective anti-ferromagnetic coupling 
term between two neighboring blocks of spins. The anti- 
phase described by \ips) is of lower energy J' than the 
in-phase by |*0t)- Notice that it is different from the 
case that two neighboring spins across the hole domain 
interacts with an anti- ferromagnetic coupling 0, H^| . 
( 8 The illustration of how to flip a block of spins is shown 
in FigEI Several mean-field type studies neglecting this 
low-energy interaction |l3| cannot distinguish anti-phase 
and in-phase. 

Because of no-double occupancy at every site, the spin 
magnitude also form a period of 2R, in which 



Now because two subspaces Hs and Ht are orthogonal 
to each other, the first-order correction to the ground 
state energy due to perturbed term Hf' a is 



H t :0t gives the transverse density p(x) in (—R,R) to the 
leading order can be obtained from the 



and repeat for a period 2R. It would be convenient to 
define qo = (<7oj0), where q n = n/(2R). The charge- 
density wave induced by stripe fluctuation is 



where Pl = 3^ and p 2 = -j^^j- 

Eas. (|10|) - (|11|) tells that two almost degenerate states 



and hence 



S(r) = \(l-p(f)) 

= S0 + S1 cos(2g -r) + S 2 cos(4g ■ r) + . . (15) 



2rt 



and So 
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where S = \{l - S),S 1 = -3^, — -2 - j^T^ 
Because of the anti-ferromagnetic interaction between 
blocks of spins across the hole stripes, we can deter- 
mine classically the z-componcnt of the spin S z (r) which 
gives 

S z {f)e^- p = Sf sm(q ■ f) + SI sin(3g ■ r) + ■ ■ ■ ( 16 ) 



_2 56 

3tt 3(2+tt) ' 

and etc vanish be- 



and S3 = 



where Q = (tt,tt), Sf = $ - 
Notice that the components S|, S*|, 
cause of the anti-phase symmetry. Including the effect 
from lattice distortion gives the same result. The non- 
vanishing higher harmonic should have wavevector 3qo, 
bqo, and etc. The schematic diagram of the classical spin 
state is shown in Fig^ 

Estimated up to the order of magnitude, the neutron 
scattering intensity is more or less proportional to 



dre- kF (S z (r)S z (0)) 



oc (5f) 2 S(k-Q + q )+S(k-Q-q ) 



+ {SI) 2 S(k-Q + 3q ) + S(k-Q- 3q ) 



(17) 



For 6<1, {S Z /SD 2 ~ 1/9. Substitute 5 = 1/8, the ratio 
is 0.093. In general, the ratio of the intensities of the sec- 
ond harmonic (Q±3qo) to that of the fundamental mode 
(Q±qo) is about 10%. However, the signal-to-noise ratio 
in current neutron scattering experiments Q is not high 
enough to observe the second harmonic peak. We expect 
further experiments on a larger pure single crystal mea- 
surement can verify our prediction. The first harmonic 
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occuring at Q±2qn vanishes. For LSCO, structural effect 
at doping concentration 5 = 1/8 enhances the fundamen- 
tal mode so that the second harmonic is even harder to 
be observed. 

In order to obtain the low- energy properties, we apply 
the linear spin- wave expansion 14], i.e., 



S z {n) = ±(S(fi) - b\h) 



(18) 
(19) 
(20) 



in which &j's and b^s describe the quantum fluctuation 
from the classical ground state. The effective Hamilto- 
nian for spin dynamics is 



H.J — En 



4 J 



(VS(n)S(n+i)bib i+1 + h.c.) (21) 

i 

In the long- wavelength limit, we replace the geometric 
mean \J S(fi)S(fi + i) by the arithmetic mean (S(ri) + 
S(fi+i))/2, and then perform the Fourier transform to 
the variables b\ and bi. Notice that the sites of classical 
spin pointing up and down are merged together in their 
Fourier modes. Then 



fc 

k 

h.c.) (22) 



Yk+2q a b k b -k-2q 



where 7^ = (cos(k x ) + cos(fc y ))/2, and the terms involv- 
ing higher harmonics are neglected. The Hamiltonian is 
quadratic and in principle can be straightforwardly diag- 
onalized. Since we are only interested in the excitation 
spectrum around k = ±2ifo, write 

HA^ ±2q - a = ]>>+ M +^+ £ J-M-V- (23) 



where r]± = (b R b^ 6 T _- bK^f and 
/ 



M± = J 



45 2S 1 iS 0% 2S l7 ^\ 

2S 1 AS Q 2S l7jE=F2? 4S 7s T2f 

V 25 l7Ma - 45 07 g T2f 25! 45 / 



(24) 



Note that it is the bosonic version of the spin-density 
wave induced by charge stripes instead of the fermionic 
one due to Fermi surface instability |15| . 



Diagonalize the Hamiltonian, we get the gapless exci- 
tation at k = ±2<fo, with the anisotropic spin- wave ve- 
locities, in which 



e fe=(fc,,0) - iS o J 



S 2 

1 ~ 7J2 sec2 (9o) 



and 



e k=(±2q ,k y ) 



4S J\k y \ 



\k x ±2q \ (25) 



(26) 



One can estimate the ratio of the anisotropic spin-wave 
velocities by the above Eas. lf^ - lt^ . 



SI 



452 sec ( qo ) 



(27) 



For LSCO, g = 2tt<5. v x /v y = 1 at 5 = 0. At 5 = 1/8, 
Vx/vy — 0.9997. Without the structural effect, one can 
safely ignore the anisotropy of spin- wave velocities. How- 
ever, perturbed structural effect can be straightforwardly 
considered by including a term 



^£cos(2(f • ri)b\bi 



(28) 



where V > enhances the spin ordering of magnitude of 
period 2R. The term behaves as a perturbed term 5M 
to enhance the first harmonic component Si of the spin 
magnitude. 



SM = 



-V 



1 

1 








1 

1 / 



(29) 



For perturbed V <C t where the weak coupling limit 
is still valid, the density profile p(x) in Ea. H12|) does not 
vary. Ecis. lf2*5l) and (|2*B1 are just modified by replacing S\ 
by Si + V/(4J). At S < 1, we can estimate 



— =M/1 



V 2 
4J 2 



(30) 



The anisotropy in the neutron scattering measurement is 
around 0.75 [l6(. Ea. H30|) gives an estimate of V ~ 1.3J. 

If the structural effect is strong enough (V/t > 1), 
for example, in the sample of the ordered stripe phase of 
La2Ni04 +< 5, we should go beyond the weak coupling limit 
such that the density profile found in Eq. (|12|l should also 
depend on V. The case of the ordered stripe phase will 
be reported elsewhere. 
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FIG. 1: Schematic diagram of spin state at doping concen- 
tration 5 = 1/8 of LSCO {q = 2n5 = vr/4). The circles 
represent the zero average of spin direction. The dashed line 
shows the envelope of the spin component in ^-direction. Any 
two neighboring "blocks of spins" feel an anti-ferromagnetic 
interaction of coupling J' . 
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FIG. 2: Illustration of the flip of a "block of spins" from (a) 
to (b). 



